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Abstract 

We discovered new hidden symmetry of the one-dimensional Hubbard model. We 
show that the one-dimensional Hubbard model on the infinite chain has the infinite- 
dimensional algebra of symmetries. This algebra is a direct sum of two s/(2)-Yangians. 
This Y(sl(2)) © Y(sl(2)) symmetry is an extension of the well-known sl(2) © sl(2) . 
The deformation parameters of the Yangians are equal up to the signs to the coupling 
constant of the Hubbard model hamiltonian. 
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1 Introduction 



Models of strongly correlated electronic systems are currently under intense study in relation 
to high-T c superconductivity. The Hubbard model is in the very center of these investigations. 
A number of results pertaining to it are now available. For a guidance in the extensive 
literature devoted to the Hubbard model we refer a reader to several reviews and reprint 
volumes [3] and the references therein. 

The first exact results in the one-dimensional Hubbard model were obtained by E. Lieb 
and F. Wu [1], who were using the Nested Bethe-Ansatz technique discovered in [2] . Since 
then, through the work of many authors, this model acquired a status of exactly solvable one. 
All of the progress achieved so far in the exactly solvable one-dimensional Hubbard model 
came through the use of the coordinate Bethe-Ansatz in combination with the si (2) © si (2) 
internal symmetry [12] exhibited by the Hubbard model hamiltonian. In particular in this 
way the exact S-matrix of the model was recently found [4]. 

Though the long-distance asymptotics of the correlaton functions in the Hubbard model 
are very well studied[5], exact formulas for them have, so far, been unaccessible by means 
of the traditional Bethe-Ansatz-based approach. In this light, in order to get the exact 
correlators and form-factors, it seems natural to turn to the symmetry-based approach which 
has been recently pioneered with impressive success in the series of works [6-8] dealing with 
the wide range of models, most notably XXZ magnet. In this approach one starts with a given 
model directly on infinite lattice and looks for an ifinite-dimensional algebra of symmetries ( 
infinite dimensional Lie algebra or some quantum group which is a deformation of an infinite 
dimensional Lie algebra) which, if large enough, permits to solve the model completely in 
the spirit close to that one of Conformal Field Theory. To carry out such a program for 
the Hubbard model one first has to identify its algebra of symmetries on the infinite chai n. 
One comes up with the result that this algebra contains a pair of commuting s/(2)-Yangians 
[9]. Each of these Yangians is an infinite dimensional quantum group which is a deformation 
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of the non-negative frequency part of the sl(2) Kac- Moody algebra, absolute values of the 
deformation parameters and the coupling constant of the Hubbard model being equal. This 
Yangian symmetry, probably, admits an extension to the Yangian Double [10], however we 
don't consider this question in the present letter. Important motivation to look for the 
Yangian symmetry in the Hubbard model came from the paper [4], where the two-particle 
S-matrix of the elementary excitations was found. This S-matrix is easily recognizable as 
an intertwiner between the tensor products of the two (simplest) evaluation representations 
[11] of the direct sum of two Yangians. This pointed out the likely existence of the Yangian 
symmetry in the model ( For the detailed general discussion of this point see, for example, 
[10] ). 

2 The Yangian symmetry of the Hubbard model 

In this section we give explicit expressions for the generators of the two Yangians commuting 
with the Hubbard hamiltonian on the infinite chain. We omit quite long and tedious ( but 
straitforward ) computations leading to the result formulated below. 

The Hubbard hamiltonian we are working with is given by the following expression: 
. _ 1 1 

i <r=T,i i 

Here c i>a are canonical fermionic annihilation operators on the lattice, % runs through the 
integer numbers and labels the sites of the infinite chain, a labels the spin degrees of freedom, 
U is the (real) coupling constant, and n ij(T = c\ a c ita is the number operator for spin a on site 
i . 

In order to formulate the result we introduce the following notations for the densities of 
the symmetry generators: 

£i = C l,1 C i+n,l i J~7 = e\,l C i+n,1 , = 4,T C *+n,T — C l,l C i+n,l ■ (2) 

Now we formulate the principal: 
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Statement 1 the six operators: (E , F , H , E±, F±, Hi) given below commute with the 
hamiltonian (1) : 

3> = E#. ^o = E-^ ffo = E«?. ( 3 ) 

i i i 

e, = E^-O-fEW-W). ( 4 ) 

i i<j 

Fi = E(^ - ^) + \ E(*W - > (5) 

i i<j 

Hx = E(^-^ _1 ) + ^E(^-^)- (6) 

Verification of this statement is a straightforward calculation. The above operators, as we 
will show in a moment, give us the generators of one of the Yangians. In order to get another 
one we proceed as follows: one easily verifies that under the canonical transformation : 

Ci,i -> c u , q iT -> (-l) l 4 T (7) 

the kinetic term in (1) doesn't change, while the interaction term changes its sign. This 
observation immediately leads to the conclusion that the operators E' , F' q , H' , E[, F[, H[ ob- 
tained from (3-6) under the transformation (7) and U — > —U commute with the hamiltonian 
as well. Thus we obtained 12 generators commuting with the Hubbard hamiltonian (1). 

Now we clarify what algebra these generators generate. First of all the operators (3) 
and their primed counterparts are the well-known sl(2) © sl(2) generators. They are known 
to commute with the Hubbard hamiltonian (under the variety of boundary conditions) on 
the finite chain also, as well as in higher dimensions [12]. The rest of the generators extend 
sl(2) © si (2) to Y(sl(2)) © Y(sl(2)) , so that all the unprimed generators commute with the 
primed ones, and between themselves the (un)primed generators satisfy the defining relations 
of the s/(2)-Yangian [9,10]. We formulate this as : 

Statement 2 i). The operators {E t , F e , H e } e=0 i provide a representation of the sl(2)- 
Yangian with the deformation parameter equal to U ( or —U ), i. e. the following defining 
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relations [9,10] are satisfied: 

[Eq, F ] = H , [H , E ] = 2E , [H , F ] = —2F . (8) 
( i.e. E ,F Q ,H generate si (2) ); 

[E ,F 1 ]=H 1 , [F ,E 1 ] = -H 1 , [H ,E 1 ]=2E 1 , (9) 

[E , Hi] = -2E 1 , [F , Hi] = 2F 1 , [H , F,] = -2F 1 , (10) 

[E ,E 1 ]=0, [F ,F 1 ]=0, [H ,H 1 ] = 0. (11) 

( Ei, Fi, Hi form a vector representation of (8) ); 
and finally one has the deformed Serre relations: 

[Hi,[Ei,Fi]] = U 2 ({H ,E ,Fi}-{H ,F ,Ei}), (12) 
[Ei,[Hi,Ei\] = U 2 ({Eo,H ,Ei}-{E ,Eo,Hi}), (13) 
[F 1 ,[H 1 ,F 1 ]] = U 2 ({F ,Ho,Fi}-{F ,Fo,Hi}), (14) 

[Hi,[Hi,Ei]]+2[Ei,[Ei,Fi]] = 

U 2 ({H , H , Ei} - {Ho, E , Hi} + 2{E , E , Fi} - 2{E , F , Ei}), (15) 
[H 1 ,[H 1 ,F 1 ]] + 2[F 1 ,[F 1 ,E 1 ]] = 

U 2 ({H , Ho, Fi} - {H , F , Hi} + 2{F , F , E x } - 2{F , E , Fj). (16) 

Where the curly brackets mean symmetrization: {xi,x 2 ,x 3 } = ^J2i 1 ^i 2 ^i 3 x ii x i 2 x i3- 

ii). The primed generators satisfy the same relations as above and commute with the 
unprimed ones. 

Now let us make some comments about the results. As one sees from the above defining 
relations, one can assign to the Yangian generated by (3-6) a deformation parameter equal 
to U or —U. To fix the sign one has to find the action of the Yangian generators upon the 
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multiparticle eigenstates of the hamiltonian i.e. the coproduct. Unfortunately in the present 
case one apparently doesn't have means to derive the form of the symmetry generators using 
the coproduct and the one-site representation of a symmetry algebra as the starting point, as 
can be done, for example, in XXZ-magnet. Hence we relied essentially upon the guesswork 
and heavy computations in order to get the symmetry generators. 

It is clear, that the expressions (3-6) provide a free-fermion representation of the Yan- 
gian. How non-trivial this representation is ? One easily sees, that taking, for example, 
commutators of the generators (E i: F 1 , Hi) beteen themselves we get operators wich con- 
tain the densities (2) with increasing absolute values of n. This apparently shows that the 
representation of the Yangian we are dealing with is infinite dimensional. Moreover, when 
U = the generators (3-6) turn into the linear combinations of the infinite-wedge (infinite- 
dimensional) [13] representation for the generators of s/(2)-Kac-Moody algebra. 

Is the above Yangian symmetry complete ? Probably it's not, and can be extended to 
the Yangian Double which is a deformation of the complete s/(2)-Kac Moody algebra [10]. 
Indeed in the case, for example, of XXZ magnet the complete algebra is a deformation of the 
complete si (2) so it is natural to look for the larger symmetry in the case of the Hubbard 
model as well. 

3 Conclusion 

As we have seen the 1-d nearest-neighbour interaction Hubbard model has the Yangian plus 
Yangian symmetry. This symmetry is an extension of si (2) © sl(2). The latter symmetry is 
generic to the Hubbard model and exists not only in the simplest case we were looking into, 
but also in higher-dimensional cases and in the Hubbard models on non-periodic lattices 
[12]. It would be interesting to see if these more comlicated Hubbard models admit some 
extensions of their si (2) © si (2) symmetries. 

Returning to the model (1) the next step is to try to utilise the Yangian symmetry in order 
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to get the correlators and the form-factors. This includes the identification of the Hilbert 
space of the model as a representation of the Yangians and expression of the elementary 
excitations' creation opertors in terms of vertex operatos related to the Yangians. 

It will be instructive to understand the relation between the Yangian symmetry and the 
abelian algebra of conservation laws discovered in [14]. 

It also would be interesting to try to modify the Hubbard model ( by the introduction 
of some non-local interaction ) in such a way that a new model would exhibit a Yangian 
symmetry on the finite lattice. This program was recently performed for the Heisenberg 
model in [15]. 
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